LETTER TO THE EDITOR 



Transverse frames for Petrov type I spacetimes: 
a general algebraic procedure 

Virginia RefJ§, Marco Brunij, David R. Matraversf and 
Frances T. Whitej 

f Institute of Cosmology and Gravitation, 
Portsmouth University, Portsmouth POl 2EG, U.K. 
\ Dipartimento di Fisica "E.R. Caianiello", 

Universita di Salerno, via S.Allende 84081 Baronissi , Salerno, Italy 
§ Istitituto Nazionale di Fisica Nucleare, Sezione di Napoli 

Abstract. We develop an algebraic procedure to rotate a general Newman- 
Penrose tetrad in a Petrov type I spacetime into a frame with Weyl scalars and 
\l/3 equal to zero, assuming that initially all the Weyl scalars are non vanishing. 
The new frame highlights the physical properties of the spacetime. In particular, 
in a Petrov Type I spacetime, setting Vti and to zero makes apparent the 
superposition of a Coulomb-type effect ^2 with transverse degrees of freedom 
and *4. 



PACS numbers: Add 
29 November 2002 



1. Introduction 

The most promising sources of gravitational waves, for example the merging of two 
black holes, are currently studied primarily using numerical methods. Using the 
standard tools of classical relativity to characterize the space-time in an invariant 
manner we can complement this numerical work, contributing to the physical 
interpretation of the output. There are many tools developed for analyzing exact 
solutions which can be applied in the context of numerical relativity. Here we will 
work within the Newman-Penrose formalism, on the problem of characterizing a Petrov 
type I spacetime. Indeed one expects that realistic gravitational wave sources will be 
described by a radiating spacetime of type I. 

In the Newman-Penrose formalism the ten independent components of the Weyl 
tensor are expressed as five complex scalars which are the projection of the Weyl 
tensor on the null tetrad {1, n, m, m}. In vacuum, all the physical information on 
the curvature is contained in these scalars. Rotations of the tetrad, conveniently 
divided into three classes, bring about transformations of these five Weyl scalars (see 
for example that can be used to produce simplifications. 

The Petrov classification of the Weyl tensor depends on the number of repeated 
roots of the equation i]/^ 6 ™' = fj, where ^i^ New ^ is obtained by a class II rotation and 
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is expressed in terms of the initial scalars and the rotation parameter Q , |Q) . When 
some of the roots coincide the spacetime is said to be of special type. A pure plane 
gravitational wave is Petrov type N (four coincident roots) and in this case the tetrad 
can be chosen to have W4 (or Wo) as the only non-zero Weyl scalar. The scalar W2 
represents a Coulomb-type tidal field. In a type D spacetime (two distinct repeated 
roots) such as Schwarzschild or Kerr a tetrad can be found such that W2 is the only 
non-zero scalar. For Petrov type II spacetimes (three distinct roots, one of which is 
repeated) a frame can be found where W4 (or Wo) and W2 are the only non-zero scalars. 
In HI a series of explicit algorithms for putting the Weyl scalars into these standard 
forms, allowing for all the possible initial non-zero components that one can have for 
each special Petrov type, is presented. 

The most general spacetime is Petrov type I, in which ^^ New ^ = has four 
distinct roots, and here the physical meaning of the Weyl scalars can be highlighted 
if one chooses the tetrad so that the scalars appear in one of the standard forms ||] : 

{0,1-1, 1-2, 1-3,0} (1) 

{W ,0,W 2 ,0,W 4 }. (2) 

The first case is the easier to achieve, at least in principle, but the second case is 
physically more interesting in the light of the physical interpretation of the Weyl 
scalars |[ outlined above. Since Wi and W3 represent longitudinal "gauge fields" 
they can conveniently be put to zero, with Wo and W4 representing transverse degrees 
of freedom, and W2 a Coulomb-type field. In the following, we call a tetrad in which 
the scalars take the form (2) a transverse frame, and we shall denote the Weyl scalars 
in this tetrad as Wj, Wj and Wj §]. 

If the spacetime is radiative the dominant Petrov type in the wave zone will be N 
(Peeling-Off Theorem M), and one expects Wj and Wj to represent gravitational 
waves. If evaluated at a finite distance from the source, WJ and Wj will still 
be transverse, but in general will include other effects. In relation to this, it is 
worth pointing out that a Petrov type I spacetime does not necessarily incorporate 
gravitational waves (for example there are stationary metrics of type I), nor does 
the presence of gravitational waves necessarily imply a non vanishing W 4 or Wo, as 
it is clear from the existence of form (1) for any type I spacetime that may include 
radiation. 

Because of the physical interest of transverse frames we develop here an algebraic 
procedure for arriving at a transverse tetrad for a generic Petrov type I spacetime. 
We give expressions for the final rotated scalars in terms of the initial scalars, which 
are all assumed to be non vanishing. This is likely to happen if a completely arbitrary 
choice of the initial tetrad frame is made. The cases in which some of the initial scalars 
are zero in Petrov type I are dealt with in Q . 

Section 2 is the main body of the letter, containing the details of the procedure. 
Section 3 consists of conclusions and discussion. The final expressions for the scalars 
that we obtain after performing the transformations are shown in the Appendix. 



2. The Procedure 



The problem for a given type I spacetime is how to arrive at a transverse frame 
{W , 0, W 2 , 0, W4} from the general case of all scalars non-zero. To do this we will use 
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class I and class II transformations, which arc jlj] 



Class I * — > *o 

*i -> *i + 0**0 

* 2 ^*2 + 2a** 1 + (a*) 2 * 

*3 -> *3 + 3a** 2 + 3(a*) 2 *i + (a*) 3 * 

* 4 -> *4 + 4a** 3 + 6(a*) 2 * 2 + 4(a*) 3 *i + (a*) 4 * 

Class II * -> * + 46*i + 66 2 * 2 + 46 3 * 3 + 6 4 * 4 
*i -»• *i + 36*2 + 36 2 * 3 + 6 3 * 4 

* 2 _> \J/ 2 + 26* 3 + 6 2 * 4 

*3 -> *3 + 6*4 
*4 -> * 4 

where a and 6 are complex parameters. 

Starting from the general situation in which all the five scalars are different from 
zero, it is not easy to obtain a transverse frame by directly carrying out the two types 
of rotations. To do this we would need to solve the resulting system of equations for 
a and 6 : 

*i T = *i + 6*o + 3a(* 2 + 26*! + 6 2 * ) 
+ 3a 2 (* 3 + 36*2 + 36 2 *i + 6 3 * ) 

+ a 3 (* 4 + 46*3 + 66 2 * 2 + 46 3 *i + 6 4 * ) = (3) 
* 3 T = * 3 + 36*2 + 36 2 *! + 6 3 * 

+ a(* 4 + 46*3 + 66 2 * 2 + 46 3 *i + 6 4 * ) = (4) 

Here a class I transformation with parameter 6 has been followed by a class 
II transformation with parameter a. The difficulty is that the degree of the above 
equations is too high for the problem to be solved by a direct approach. In order 
to get around this, we perform preliminary rotations on the original scalars to put 
*4 = and *o = 0, arriving at the standard form (I). This simplifies the system 
considerably because *o is the coefficient of 6 4 and *4 is the coefficient of a 4 in the 
transformation equation for the *'s. So these preliminary rotations ensure that the 
final equations to be solved are only third order. 

In total we perform four rotations: 

^ Class I ~j Class II ^ ff Class I Class II \J/~T 

where we require: 

* = {*o,*i,* 2 ,*3,*4}, 

*' ={* ,*' 1 ,*2,*3,0}, (5) 

*" ={0,*i',* 2 ',*3',0}, 
* T = {* T ,0,* 2 r ,0,*J}. 

In order to get * 4 = 0, we apply a class I rotation, with parameter 6, to the 
original scalars. The expression for * 4 is a fourth order polynomial involving all the 
original scalars. This can be solved by Maple. Substituting one of the solutions for 
6 into the rotation equations, we obtain expressions for the rotated scalars in terms 
of the original ones. These turn out to be extremely long and must be manipulated 
before they can be used in the further calculations. 
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The next step is to set \l/g equal to zero. We apply a class II rotation to the 
scalars (under which \& 4 ' = ty' 4 remains zero) and require = 0. This means 
solving a third order polynomial equation in a. Choosing one of the solutions for a, 
we get the resulting non-zero scalars ('I'", ^2, ^3) in terms of the old ones. 

Finally we rotate twice more, first class I (with parameter g) and then class II 
(with parameter /). Since in our starting frame \&q and \&4 are zero, the resulting 
equations for the new iffj and tyj are as follows: 

* X T = + 3/(*2 + 2p#i) + 3/ 2 (* 3 ' + 30*2 + 3g 2 V'{) 

+ f(4g^ + 6gH' 2 ' + 4g 3 H>'() = (6) 

* 3 T = *3 + 3.g*2 + 3.g 2 *i + f{±9% + 6.g 2 *i' + 4g 3 *") = (7) 
In this frame, the equations are only third order in both parameters and the system 
is straightforward to solve. The solution is formally simple. For example the possible 
solutions include || : 

9 = \l^r» f = 1 /9 (8) 

However, the full expressions of the above quantities, in terms of the original 
scalars "J, are too long even to be displayed. To make them more manageable, 
we implement the Maple command optimize. This telescopes the expressions by 
assigning to each sub-expression a name: ti,t%,..., etc and builds up a chain of 
definitions relating these new objects. This greatly simplifies the implementation 
of these expressions in computer language. 

The final expressions for the scalars and all the definitions required for their 
compact expression are reported in the Appendix. Here, for the sake of brevity, we 
report only the last line for each scalar: 

^0 = — 2ti4o/£i43 + 3/2^132^140 ($2 + *91 + *92 + *149 + 2*143*140) 

- 1/2(12*79 + 12* 82 + 4* 86 + 3*i 43 *i6o)/*i55 

+ 1/16*167(4*143*87 + 6ti4 2 ii60 + 4*i55*i4o)/*97 

(9) 

= *2 + *80 + *82 + *136 + 2*145*142 

— (12*87 + 12*89 + 4*93 + 3*145*152)/*145 

+ 1/4*132*142(4*145*94 + 6*144*152 + 4*145*144*142) 

(10) 

^4 — 4*143*87 + 6*142 (^ + *91 + *92 + 2*134*37) + 4*i43*i 42 *i40 

(11) 

where the t's are temporary objects created by Maple, built from the original scalars. 
3. Conclusions 

The transverse frame that we arrive at after our four rotations is not unique, as there 
is a two-fold infinity of tetrads satisfying ^4=^/3=0. This arbitrariness corresponds to 
the extra freedom that one has of performing a further rotation of class III, depending 
on two real parameters, in order to completely fix the tetrad ||. 

The expressions for ^2, ^0 an d ^4 we obtain are quite complicated (see the 
Appendix) . Partly they depend on the specific root of the polynomial one chooses at 
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each step of the procedure (|5|) . We have not investigated here wheater one particular 
path in this general scheme may bring about more manageable expressions for the 
scalars. In any practical case one will have to substitute the expression for the initial 
Weyl scalars in terms of the given metric and its derivative and these will be the 
expressions that one will want to simplify. 

Nevertheless our results are completely general and can be used to put any Petrov 
type I spacetime into the transverse frame in the case when the initial scalars are all 
non vanishing. This is useful whenever one wants to investigate the physical properties 
of a Petrov Type I spacetime using curvature variables. A good example is the full 
non-linear numerical treatment of a perturbed black hole. In this case the transverse 
frame is the best way to highlight the perturbations affecting the curvature held, which 
in the unperturbed case is represented by x &2 only. 

More generally, it is of interest and potential value that our optimized expressions 
can in principle be applied to any numerically computed spacetime in which one 
produces numerical values for the initial scalars. 
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4. Appendix 

Final expression for the rotated scalar ^/q T : 

7(6), 
*?, 

*4 * #2, 
£8 * * 2 . 

v/(3), 

* 4 *tl3, 

£8 2 , 
t2 2 , 

y/l^i * {-tl2 * £13 * *o + 18 * £16 * t8 - 81 * * * tl9 - 54 * £3 * *o * £2 + 
27**4* £25 + 54* £2* £9)), 

(-27 * £3 * * + 27 * * 4 * t2 + 27 * £9 + 3 * til * £32) 1/3 , 
t2 * £36, 



£1 = 
£2 = 
£3 = 
£8 = 
£9 = 
£11 = 
£12 = 
£13 = 
£16 = 
£19 = 
£25 = 
£32 = 
+ 

£36 = 
£37 = 
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£39 = * 2 * *o, 
£40 = £39 * £36, 
£42 = £36 2 , 
£48 = l/£36, 

£50 = VT(6 * £37 - 6 * £40 + £42 * * + 3 * £16 + 9 * * * *8) * £48), 
£77 = VT(72 * £37 * £50 - 72 * £39 * £36 * £50 - 6 * £50 * £42 * * - 

+ 18 * £50 * * 4 * £13 - 54 * £50 * * * t8 - 72 * £2 * *i * £1 * £36 + 

+ 108 * * x * £1 * £40) * £48/£50), 
£78 = £1 * £50 + £77, 
£79 = 1/6* £78**2, 
£81 = 1/36 *£78 2 , 
£82 = £81 * *i, 
£86 = 1/6* £81* £78**o, 
£87 = 3 * £79 + 3 * £82 + £86, 
£88 = 1/6* £78**o, 
£89 = *i +£88, 
£90 = 1/6* £78**i, 
£91 = 2* £90, 
£92 = £81 * * , 
£93 = * 2 + £91 + £92, 
£94 = £89 * £93, 
£97 = £87 2 , 

£100 = £93 2 , 
£101 = £100* £93, 
£103 = £89 2 , 

£117 = VT64 * £103 * £89 * £87 - 36 * £103 * £100 - 108 * £94 * £87 * * + 

+ 27 * £13 * £97 + 54 * * * £101), 
£122 = (54 * £94 * £87 - 27 * * * £97 - 27 * £101 + 3 * £11 * £117 * £87) ( l/3), 
£123 = 1/£122, 
£132 = l/£87, 

£134 = - 2 * £89 * £123 + (1/6 * £122 + 3/2 * £100 * £123 - 1/2 * * 2 - £90 - 1/2 * £92) * £132, 
£137 = £134 2 , 

£140 = *i + £88 + 3 * £134 * £93 + 3 * £137 * £87, 

£142 = £87/£140, 

£143 = \f{t!42), 

£149 = 2 * £134 * £87, 

£155 = £143* £142, 

£160 = * 2 + £91 + £92 + £149, 

£167 = £140 2 , 
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£178 = - 2 * £140/£143 + 3/2 * £132 * £140 * (* 2 + *91 + £92 + £149 + 2 * £143 * £140) - 
+ 1/2 * (12 * £79 + 12 * £82 + 4 * £86 + 3 * £143 * £160)/£155 
+ 1/16 * £167 * (4 * £143 * t87 + 6 * £142 * £160 + 4 * £155 * £140)/£97 

Optimized expression for the rotated scalar * 2 T : 

tl - 7(6), 
£2 = * 2 , 

£3 = * 4 ** 2 , 

£8 = * 2 , 

m = £8** 2 , 

til =7(3), 

£12 =*2, 

£13 = * 2 , 

tl6 = * 4 *£13, 

£19 = £8 2 , 

£25 = £2 2 , 

£32 = v/(*4 * (— tl2 * £13 * * + 18 * £16 * £8 - 81 * * * tl9 + 

- 54 * £3 * * * t2 + 27 * * 4 * £25 + 54 * £2 * £9)), 

£36 = (-27 * £3 * * + 27 * * 4 * £2 + 27 * t9 + 3 * til * £32) 1/3 , 

£37 = £2* £36, 

i39 = * 2 * *o, 

MO = £39* £36, 

£42 = £36 2 , 

£48 = l/£36, 

£50 = VT( 6 * £37 - 6 * £40 + £42 * * + 3 * £16 + 9 * * * *8) * £48), 

£77 = \f{{72 * £37 * £50 - 72 * £39 * £36 * £50 - 6 * £50 * £42 * * + 

- 18 * £50 * * 4 * £13 - 54 * £50 * * * t8 - 72 * £2 * *i * £1 * £36 + 





+ 108 * *i * £1 * £40) * £48/£50) 


£78 


= £1 * £50 + £77, 


£79 


= 1/6* £78**i, 


£80 


= 2 * £79, 


£81 


= 1/36 *£78 2 , 


£82 


= £81 * * , 


£83 


= 1/6* £78**o, 


£84 


= #i +£83, 


£85 


= * 2 + £80 + £82, 


£86 


= £84 * £85, 


£87 


= 1/6* £78** 2 , 
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tS9 =£81**1, 

£93 = 1/6 * £81 * £78 * * , 

£94 = 3 * £87 + 3 * £89 + £93, 

£97 = £94 2 , 

£100 = £85 2 , 

£101 = £100* £85, 

£103 = £84 2 , 

£117 = VT64 * £103 * £84 * £94 - 36 * £103 * £100 - 108 * £86 * £94 * * + 

+ 27 * £13 * £97 + 54 * * * £101), 
£122 = (54 * £86 * £94 - 27 * * * £97 - 27 * £101 + 3 * £11 * £117 * £94) 1/3 , 
£123 = 1/£122, 
£132 = l/£94, 

£134 = - 2 * £84 * £123 + (1/6 * £122 + 3/2 * £100 * £123 - 1/2 * * 2 - £79 - 1/2 * £82) * £132, 
£136 = 2 * £134 * £94, 
£139 = £134 2 , 

£142 = #i + £83 + 3 * £134 * £85 + 3 * £139 * £94, 

£144 = £94/£142, 

£145 = VT*144), 

£152 = * 2 + t80 + £82 + £136, 

£168 = * 2 + £80 + £82 + £136 + 2 * £145 * £142 - 

+ (12 * £87 + 12 * £89 + 4 * £93 + 3 * £145 * £152)/£145 + 

+ 1/4 * £132 * £142 * (4 * £145 * £94 + 6 * £144 * £152 + 4 * £145 * £144 * £142) 

Optimized expression for the rotated scalar ^4 : 

£1 = v/(6), 
£2 = * 2 , 

£3 =*4** 2 , 

£8 = 

£9 =£8**2, 

£11 =vf3), 

£12 =*2, 

£13 = * 2 , 

£16 =* 4 *£13, 

£19 =£8 2 , 
£25 =£2 2 , 

£32 = y/l^i * (-£12 * £13 * * + 18 * £16 * £8 - 81 * * * £19 + 

- 54 * £3 * * * t2 + 27 * * 4 * £25 + 54 * £2 * £9)), 
£36 = (-27 * £3 * * + 27 * * 4 * £2 + 27 * £9 + 3 * £11 * £32) 1/3 , 
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£37 = £2* £36, 
£39 = * 2 * *o, 
£40 =£39* £36, 
£42 = £36 2 , 
£48 = l/£36, 

£50 = VT(6 * £37 - 6 * £40 + £42 * * + 3 * £16 + 9 * * * *8) * £48), 
£77 = VT(72 * £37 * £50 - 72 * £39 * £36 * £50 - 6 * £50 * £42 * * + 

- 18 * £50 * * 4 * £13 - 54 * £50 * * * *8 - 72 * £2 * * a * £1 * £36 + 

+ 108 * *i * £1 * £40) * £48/£50), 
£78 = £l*£50 + £77, 
£79 =1/6* £78**i, 
£80 = 2 * £79, 
£81 =l/36*£78 2 , 
£82 = £81 * * , 
£83 =1/6* £78**o, 
£84 = *i + £83, 
£85 = * 2 + £80 + £82, 
£86 =£84* £85, 
£87 =1/6* £78**2, 
£89 =£81**i, 
£93 = 1/6* £81* £78**o, 
£94 = 3 * £87 + 3 * £89 + £93, 
£97 =£94 2 , 
£100 = £85 2 , 
£101 = £100* £85, 

£103 = £84 2 , 

£117 = v/(64 * £103 * £84 * £94 - 36 * £103 * £100 - 108 * £86 * £94 * * + 

+ 27 * £13 * £97 + 54 * * * £101), 
£122 = (54 * £86 * £94 - 27 * * * £97 - 27 * £101 + 3 * £11 * £117 * £94) 1/3 , 
£123 = 1/£122, 
£132 = l/£94, 

£134 = - 2 * £84* £123+ (1/6 * £122 + 3/2 * £100 * £123 - 1/2* * 2 - £79 - 1/2 * £82) * £132, 
£136 = 2* £134* £94, 
£139 = £134 2 , 

£142 = *i + £83 + 3 * £134 * £85 + 3 * £139 * £94, 

£144 = £94/£142, 

£145 = \f{t\4A), 

£152 = * 2 + £80 + £82 + £136, 

£168 = * 2 + £80 + £82 + £136 + 2 * £145 * £142 - 
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+ (12 * £87 + 12 * £89 + 4 * £93 + 3 * £145 * £152)/£145 + 

+ 1/4 * tl32 * il42 * (4 * £145 * £94 + 6 * £144 * £152 + 4 * £145 * £144 * £142) 
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